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EFFECT OF AXIAL LOAD ON MODE SHAPES AND FREQUENCIES OF BEAMS 

by Francis J. Shaker 
Lewis Research Center 


SUMMARY 

An investigation of the effect of axial load on the natural frequencies and mode 
shapes of uniform beams with various types of boxmdary conditions and of a cantilevered 
beam with a concentrated mass at the tip is presented. This investigation yielded ex- 
pressions for the mode shapes and characteristic equations for the various cases con- 
sidered. For the uniform beams the characteristic equations were solved either numer- 
ically or in closed form, and the results are presented by a series of graphs showing the 
effect of preload on the different types of beams. The effect of axial load on the mode 
shapes is also shown in graphical form for several different loading conditions. For the 
cantilevered beam with a tip mass, two types of axial loads were considered: In the 
first case the axial load vector remained constant, and in the second case the load was 
directed through the root of the beam at all times. The results of this portion of the in- 
vestigation are presented in graphs that show the effects of both tip -mass variation and 
axial load on the fundamental frequency of the system. 


INTRODUCTION 

In the design of certain spacecraft structural components it sometimes becomes 
necessary to determine the normal modes and frequencies of beam t 3 rpe components 
which are in a state of preload or prestress. This preload may result from inertial ef- 
fects, as in the case of a spinnii^ spacecraft, or it may be induced by mechanical 
means. For example, current designs for large, flexible solar arrays (refs. 1 to 4) 
are such that the boom that supports the array is in a state of prestress due to the ten- 
sion that must be maintained in the solar cell substrate. Also, certain types of attach- 
ment and ejection mechanisms can cause preload in parts of the primary spacecraft 
structure. In most of these cases an exact solution for the normal modes and frequen- 
cies of the system would be intractable and recourse would be made to one of the many 
approximate solutions such as the Rayleigh-Ritz method, Galerkin’s method, finite 
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elements, etc. (refs. 5 to 8). There exists, however, a class of problems for which an 
exact solution can be obtained. This class consists of xmiform beams under constant 
axial load with various types of simple boundary conditions. This t^pe of structure may 
represent some structural components reasonably well for preliminary design informa- 
tion. Even in those cases where it does not adequately model the structure, the modes 
obtained for these cases can be used to obtain Rayleigh -Ritz or Galerkin type of 
solutions. 

For these reasons an investigation was made to determine the effect of an axial load 
on the modes and frequencies of a uniform beam with various boundary conditions and of 
a cantilevered beam with a tip mass. The boundary conditions considered for the uni- 
form beam include all possible simple boundary conditions at both ends of the beam 
(i.e. , free, hir^ed, clamped, and guided). For the cantilevered beam with a tip mass, 
two types of axial loads were considered. The first type was a constant axial load ap- 
plied at the tip, and the second was an axial load whose direction is always throv^h the 
fixed support. The second is the limiting case for some current large solar array de- 
signs if the mass of the blanket is negligible compared with the mass of the support 
beam. When the mass of the blanket is not negligible, which is usually the situation, 
these modes could still be used in a Rayleigh-Ritz analysis for this type of array. 


SYMBOLS 


A,B,C,D,E 

■^n’^n 

El 

k 

k 

cr 


I 

M(x,t),m(x), 

m(x) 


arbitrary constants of integration 

•Mn 

constants of integration corresponding to mode of vibration 
\aniform beam bending stiffness 
axial load parameter, |/p/EI 
nondimensional axial load parameter, yPz/EI 

no ndimensio nal axial load parameter corresponding to critical load, 
length of beam 

moment distribution along length of beam 


Jl-Q 



total mass of beam, pi 

tip mass for cantilevered beam 

axial load 

critical or Euler ’s buckling load for beam 
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Q(x,t),q(x), 

q(x) 

t 

V(x,t),v(x) 

v(x) 

X 

X 

“2n 


Vrj, 

P 

(V 

Subscripts: 


shearing force distribution 
time 

beam displacement 

nondimensional beam displacement, v/l 
th 

mode shape of n‘'" mode 
lengthwise coordinate 
nondimensional lengthwise coordinate, x/l 
characteristic values, eq. (9) 
characteristic values for mode 
frequency parameter, ypw/EI 
nondimensional frequency parameter, 01 

nondimensional frequency parameter corresponding to n*^" mode of 
vibration , 

ratio of tip mass to total beam mass, 
mass per unit length of beam 
circular frequency of vibration 
circular natural frequency of n*^" mode 


B 

n 

T 


beam 

integral number designatii^ natural bending mode and frequency of beam 
tip of cantilevered beam 


Superscript: 

differentiation with respect to x 


THEORETICAL ANALYSIS 
Uniform Beam Under Axial Load 

Equ ations of motion. - The development of the equations governing the bending 
vibrations of uniform beams are presented in numerous texts on vibration theory 
(refs. 5, 6, or 8). In the development that follows the effect of an axial load is included. 
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Consider a uniform beam vibrating freely vmder the action of a constant axial load 
P as shown in figure 1(a). Let V(x, t) be the displacement of any point x along the 
neutral axis of the beam. The internal shear, moment, and inertial forces acting on an 
element dx in length will be as shown in figure 1(b). Summing forces and moments on 
this element yields the following two equilibrium equations; 

3Q(x,t) ^ 3^V(x,t) . „ 

3x g|.2 

+ P 9V(M) ^ ^ Q 

ax ax 

From elementary beam theory the moment-curvature relationship for the sign convention 
depicted In figures 1 is 


M(x,t) = El 


aV(x,t) 

ax2 


(10 


The shearing force Q(x,t) can be expressed in terms of the displacement variable by us- 
ing equations (lb) and (Ic). Thus 


Q(x, t) 


El + P 5V(x,t) 

ax^ 


( 2 ) 


Finally, from equations (la) and (2) the equation of motion in terms of the displacement 
variable can be written as 



P a^V(x,t) 


El 


ax'^ 


p dh(x,t) 

El 3,2 


= 0 


(3) 


Equation (3) represents the governing partial differential equation describing the motion 
of the beam, and equations (Ic) and (2) represent the moment and shear distribution 
along the beam. 

Solution of equation and boundary conditions . - The time variable in equations (1) to 
(3) can be eliminated by assuming a solution of the form 

V(x,t) = v(x) sin o)t (4a) 
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M(x,t) = m(x) sin cot 


(4b) 


Q(x,t) = q(x) sin wt 


Substituting equations (4) into equations (1) to (3) yields 

= 0 

dx*^ dx^ 


m(x) = El 

dx2 


q(x) = -El 


d^v(x) ^ jj2 dv(x) 


dx*" 


dx 


(4c) 


(5) 

(6a) 


(6b) 


where = P/EI and i3^ = pch^/EI. It will be convenient at this point to put equa- 
tion (5) into nondimensional form by letting x = x/l, v = v/l, k = ki, and /3 = /3Z. Us- 
ing these expressions allows equations (5) and (6) to become 

v”"(x) + k^v**^) - F^v(x) =0 (7) 


m(x) = 


El-ti 

— V 

I 


q® = - 




( 8 ) 


The solution to equation (7), which can be readily verified by direct substitution, is 
given by 


v(x) = A cosh a jX + B sinh otjX + D cos 0 ! 2 X + E sin of 2 X 


where A, B, C, and 


D are constants and 


(9) 
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1/2 


a 


1 ~ 



"2 ■ 








(10) 


The frequency equation can now be determined from equation (9) once the boundary con- 
ditions are stipulated. The boundary conditions being considered at either end of the 
beam are as follows: 

Pinned - 


Clamped - 


V = 0 and m = 0 or v” = 0 


(11a) 


Free - 


V = 0 and v* = 0 


(11b) 


m = 0 or v’* = 0 


and q = 0 or 


v”’ + = 0 


(11c) 


Guided - 


v’ = 0 and q = 0 or v'” = 0 (lid) 

The various combinations of these boundary conditions at both ends of the beam are 
listed in table I, If the boundary conditions given by equations (11) are used at x = 0, 
two of the unknown constants in equation (9) can be evaluated, and the resulting solutions 
are for the beam free at x = 0 


r . “1 


:r . "2 


v(x) = A |cosh QijX + — cos 0^2^ )■'■ ® [sinh ckjX + — sin Q! 2 X j 

4 


O', 


(12 a) 


for the beam pinned at x = 0 


v(x) = A sinh a ^^x + E sin oigX 


(12b) 
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for the beam clamped at x = 0 


_ I ^ ca^ 

v(x) = A (cosh Q!jX - cos o-gx) + B I sinh a^x - _i sin a; 2 X 


Ofr 


and for the beam guided at x = 0 

v(x) = A cosh oijX + D cos ofgX 


(12c) 


(12d) 


Note that each of equations (12) contains two unknowns. When the boundary conditions at 
X = 1 are specified, these equations will yield two homogeneous, algebraic equations in 
these two unknowns. The determinant of the coefficients of the unknowns must be equal 
to zero for a nontrivial solution to exist. Equating this determinant to zero then yields 
the characteristic equation from which the frequency parameter /3 can be determined as 
a function of axial load P or axial load parameter k. To illustrate, consider the case 
where the beam is free at jc = 1. From equation (11c) the bovindary conditions are 

v"(l) = 0 v”’(l) +kV(l) = 0 (13) 


From equations (12c) and (13) the followii^ two homogeneous, algebraic equations in A 
and B are obtained after some manipulation; 

A^a^ cosh - Q!j cos " o'! “2^ “ ® 

A^ofg sinh + 0 !j sin 0 ^ 2 ^ ®1 " *^2 “l^~ 

Equating the determinant of the coefficients of A and B in equations (14) to zero, ex- 
panding, and simplifying will then sdeld the following characteristic equation for the free- 
free beam: 


(14a) 

(14b) 


2/3®(l - cosh CKj cos ffg) +ic^^ + 3/3^) smh sin CX2 = 0 (15a) 

Similarly, for the other boundary conditions at x = 1 the following characteristic equa- 
tions are obtained: for the free-guided beam 


3 3 

0*2 sinh otj cos 012 + cxj cosh sin oi^ 


= 0 


(15b) 
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L 


for the free -pinned beam 


3 3 

-a2 cosh sin 0 f 2 + sinh cos «2 “ ® 
for the clamped -free beam 

2j^ + (2^^ + k^) cosh cos Q !2 - sinh oij sin 012 ~ ^ 
for the clamped -pinned beam 

Q-j cosh Q!j sin Q !2 - a2 sinh cos 02 = 0 
for the clamped -guided beam 

Q!j sinh ofj cos a2 + 0(2 cosh sin <^2 = 0 
for the clamped -clamped beam 

o — 2 

2)3 (1 - cosh Q!j cos q! 2 ) “ ^ sinh otj sin «2 ~ ® 
for the guided -guided and pinned -pinned beams 

sin ^2 = 0 

and for the guided -pinned beam 


(15c) 


(15d) 


(15e) 


(15f) 


(15g) 


(15h) 


cos ^2 “ ® 

Equations (15g) and (15i) yield the following expressions for the frequency parameter: 
for the pinned-pinned or guided-guided beam - 


)3^ = mr V(ns-)^ - k® 

and for the pinned -guided beam 



(16) 


( 17 ) 
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Equations (16) and (17) exhibit the main characteristics of alL the frequency equations 
given by equations (15). Namely, as P (or k) increase, the frequencies decrease tmtil 
P eqTials the Euler buckling load V of the beam. At this point the lowest elastic fre- 
quency becomes zero. For P greater than P , the fundamental frequency becomes 
complex, and the corresponding mode shape is unstable (ref. 6, p. 451 and ref. 8, p. 
302). Thus, the maximum value for P of practical interest is P^^. The values of 
P and k for beams with the different botmdary conditions are tabulated in table II. 

wX V/X 

To determine the frequency parameter for the beams with other boundary conditions, 
equations (15) can be solved numerically for for any value of k less than the cor- 
responding kg^. 


Cantilevered Beam Under Axial Load With Tip Mass 

For the case of a cantilevered beam with a concentrated mass at the free end, two 
types of axial loads will be considered, hi the first case (fig. 2(a)) the axial load will be 
taken as a constant in both magnitude and direction so that this case corresponds to the 
uniform clamped -free beam case except that a concentrated mass is added to the free 
end. Ih the second case (fig. 3(a)) the axial load is assumed to be directed throi^h the 
root or fixed end of the beam at all times. For this case the axial load is not constant 
since it charges direction during the motion. This type of loading will be referred to as 
a directed axial load. For both cases the deflection is given by (eq. (12c) since the 
boundary condition at the fixed end is the same as for a uniform cantilevered beam. 

Also, tiie moment boundary condition at the free end is the same for both cases and is 
given by 


P’(l) =0 (18) 

The difference between the two eases will be in the shear boundary condition at the free 
end. For the ease of constant axial load, the shear botmdary conditions can be found by 
satis^ng the equilibrium condition at the free end which is depicted in figure 2(b). 

Thus 

q(l) + J?^o)^v(l) = 0 (19) 

From equations (6) and (9) this boundary condition becomes 

v'”(l) -i- k^F(l) -f 7j,j,^v(l) = Q (20) 
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v^here rjrj, is the ratio of the tip mass jKrj, to total beam mass, = pS.. From equa- 
tions (12c), (18), and (20) the following two homogeneous equations in A and B are 
obtained: 


>( 21 ) 


cosh 0 !j + Og cos sinh ckj sin q! 2 ^ = 0 

a|ck 2 sinh sin or^ + ? 7 ,pi 0 ^(cosh - cos OTg)!^ 

+ b|q!^ cosh Q!j + Q!j COS + 77r£.^(0!2 ^ 


Equating the determinant of the coefficients of A and B to zero, expanding and simplify- 
ing yields the following characteristic equation for a cantilever beam with a tip mass 
under constant axial load. 

-2jS^ + sinh sin - (2^^ + k^) cosh cos 


^T^(“l “2)^*^! ®1 *^2 “ “2 ®1 ®2^ “ ® 


( 22 ) 


For the case of a cantilevered beam with a tip mass and a directed axial load, the con- 
dition at the free end is shown in figure 3(b). From this figure the equilibrium condition 
can be written as 


v"’(l) + k^v'(l) + (vrj,^ - k^^v(l) = 0 (23) 

—4 —2 

Comparing this equation with equation (20) shows that the term 7jrp/3 - k replaces 

—4 —2 —4 ^ —4 —2 

q,p/3 . Thus multiplying equation (22) by /3 and replacing 7]rp/3 with (Tjrp/3 - k ) 3 ^elds 

the characteristic equation for the beam with a directed axial load. 

-2^ + /3 k^ sinh sin 02 - (2/3 + k ) cosh cos 012 

+ “1 “2 " ^2 ^1 °'2^ ~ ® 

Mode Shapes For Beams Under Preload 

The mode shapes for abeam under preload can be determined from equations (12) 
and one other boundary condition. For example, assume that the n^^ root of the fre- 
quency equation is known and is given by Correspondii^ to there will be an 
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and which are obtained by substitutii^ /3_ into equations (10). For the case 

in ZTi fu n 

of a free -free beam, the displacement in the n mode is given by equation (12a) as 


^n<^) = ^n 




cosh QtjjjX - 


r In 


cos 


“2n^ 


+ B, 


n 


I sinh + 


. ®2n 


a, 


a 


2n 


In 


sm 


(25) 


For the free-free beam case either of equations (14) can be used to determine a relation 
between and Using equation (14a) yields 




-A. 


(cosh ajjj - cos agjj) 


n 


(26) 


sinh a 


a, 


2n 


In 


sin a„ 


a 


In 


Substituting equations (26) into (25) and arbitrarily setting = 1 will then 3 deld the 
desired mode shape or eigenfimction for the free-free beam. In a similar manner the 
mode shapes for the other boundary conditions can be determined. These functions are 
summarized in the followii^ equations: for free-free and free-pinned beams 


— — I — ® In 

Vn(x) = |cosh Q!jjjX + cos Q!2n^ ' 


a, 


2n 


cosh - cos 0 f 2 n 


sinh a 


V 


In 


3 

«2n 


sm a 


2n 


a 


In 


y 


[sinh O'- X + sin x ' 

^ ! 


(27a) 
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for a free-guided beam 


— - - ® M — 

v^(x) = cosh cos 

“L 


/sinha^+^^slnorjA 


a. 


2n 


cosh ffjjj - cos 


I - ^9 - 

(sinh Of j^x + — ^ sin «2n^ 

V "in > 


for a clamped-free beam 


Vn(x) = (cosh 0 !^„x - cos a^^x) 


In^ 


"2n" 


/ Qlo \ 

/cosh cos « 2 „\ 

■4 


V 


“2n 

Sinh o!- + -ffi sin 

“in y 


/ _ a, 

sinh CKi„x iii sin a«_x 

V "2n ^ / 


for clamped -pin and clamped -clamped beams 


Vn(x) = cosh Qfjj^x - cos C!2n^ 


Lsinh a 


In 


sin a 


a, 


2n 



cosh cos 

A “2n 


for a clamped -guide beam 


Vn(x) = cosh Q!-,„x - cos oi^^x - 


In" 


2n" 


i2 1 [sinh tXj^x ri sin aoj,^ 

cosh o/j^- cos «2n /\ «2n 


for a guided^uided beam 


Vjj® = cos (nffx) 


(27b) 


(27c) 


(27d) 


I (27e) 


(27f) 
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for a guided -pinned beam 


= cos 




and for a pinned -pinned beam 


Vn(x) = sin (nnx) 


(27g) 


(27h) 


It can readily be shown that the mode shapes for the cantilevered beams with a tip mass 
are the same as those given by equation (27c). 


RESULTS AND DISCUSSION 
Uniform Beam Under Axial Load 

To determine the effect of axial load on the natural frequencies of a imiform beam, 
the characteristic equations (given by eqs. (15a) to (15g) were solved numerically by the 
method of bisection (ref. 9). In this investigation the axial load ratio was varied 

between -1 and 1, and the first three lowest frequencies were determined. The results 
for the fundamental elastic body frequencies are presented in figure 4. The pinned- 
pinned, guided -guided, and pinned -guided curves in this figure were developed from 
equations (16) and (17). In this figure the curves for the free -free and clamped -clamped 
beams appear as a single curve. Although the frequencies for these two cases were not 
identically equal over the entire range of P/Pgj.> the differences could not be shown for 
the scale used in the figure. The same is true for the free-guided and clamped -guided 
cases. The variation in frequency as a fvmction of axial load for the first few lower fre- 
quencies for each uniform beam case is given in figure 5. It can be noted that the var- 
iation in the second and third frequency for all cases shown is almost linear for the 
range of P/P_^ considered. These higher mode frequencies will go to zero when the 
second and third buckling loads are reached. However, these portions of the curves are 
of no practical interest unless the beam is prevented from buckling in its fundamental 
mode. 

To show the effect of axial loading on the mode shapes of the beam, equations (27) 
were used to determine the modes for P/P__ = -1. 0, 0, and 0. 8. These calculations 
were normalized, and the residts are shown in figures 6 to 12 for the first three modes. 
These figures indicate that the effect of axial load is greatest on the fundamental mode 
and that the effect decreases rapidly as the mode number increases. The most 
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prono\anced effect occurs for the free -pinned beam in its fimdamental mode (fig. 8 (a)). 
Note, however, that all the cases exhibit this effect as approaches 1. 


Cantilevered Beam With Tip Mass 

For a cantilevered beam with a tip mass equation (22) was solved for the case of a 

constant axial load and equation (24) for an axial load directed through the root of the 

beam. The method of solution of these equations was the same as that used for the vini- 

form beams. The effects of both tip mass ratio, rjrj^ = and axial load ratio 

P/P „ on the fimdamental frequencies were determined, and the results are shown in 

cr _2 I 

figures 13 to 16. Figure 13 shows the variation in the frequency parameters jSj yl + 77 ^ 

as a function of tip mass ratio, 77 ^,, for values of P/P^j. between -1 and 0. 8 . Note that 

^ + Tjj, equals ojj |/(ut'-j. + so that the frequency parameter for this case 

is related to the total mass of the system and not just the beam mass. Figure 14 shows 
the same data plotted as a function of axial load ratio with the tip mass ratio as a param- 
eter, Figures 15 and 16 are a similar set of graphs for the cantilever with a directed 
axial load. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
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TABLE I. - BOUNDARY CONDITIONS FOR UNIFORM BEAM UNDER CONSTANT AXIAL LOAD 


Free-free 

P 1 

X 


p 

V"(0) = 0 

v”(l) = 0 


0 

1 


v'”(0) + kV’(O) = 0 

v"’(l) + k^v' (1) = 0 

Free -guided 

P 



% p 

v"(0) = 0 

v'(l) = 0 





=1 

v’"(0) + k^’(O) = 0 

v'”(l) + kV(I) = 0 

Free -pinned 

P 

— -X 



v"(0) = 0 

v(l) = 0 






v”’(0) + k^'(O) = 0 

v"(l) = 0 

Guided- guided 




v’(0) = 0 

v’(l) = 0 


vC 



V'” + k^v'(O) = 0 

V"' + k^v'(l) = 0 

Guided- pinned 

P 5 

N 


p 

?(0) = 0 

v(l) = 0 


- 


v"'(0) + k^v'(O) = 0 

v”(I) = 0 

Clamped- free 

P 



p 

v(0) = 0 

v”(l) = 0 






v'(0) = 0 

v'”(l) + k2v'(l) = 0 

Pinned- pinned 

P 




v(0) = 0 

v(l) = 0 





v”(0) = 0 

v”(l) = 0 

Clamped- pinned 

P 1 

-^x 


p 

v(0) = 0 

v(l) = 0 




tMtt 

v’(0) = 0 

v"(l) = 0 

Clamped-guided 

l_ 

P 1 

— ^ X 

c 

p 

v(0) = 0 

v'(l) = 0 


- 

UJI 

0 


v’(0) = 0 

v'”(l) +"kV(l) = 0 

Clamped- clamped 

1 

— X 


p 

v(0) = 0 

7 ( 1 ) = 0 



7 - 5 -^ 

1 


v'(0) = 0 

v’(l) = 0 


firm, 

mr 
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TABLE II. - BUCKLING CHARACTERISTICS OF SIMPLE BEAMS 


I 


Boundary conditions 

Characteristic 

equation 

Eigenvalue, 

Free-free 

sin = 0 

Tin 

Free -guided 

cos k =0 
n 

(2n - 1) „ 
2 

Free- pinned 

sin = 0 

Tin 

Guided- guided 

sin = 0 

Tin 

Guided- pinned 

cos k^ = 0 

(2n - 1) , 
2 

Clamped-free 

cos k =0 
n 

(2n - 1) „ 
2 

Pinned- pinned 

' sin k^ = 0 
n 

njT 

Clamped- pinned 

' tan k^ = k^ 
n n 

4.493 

Clamped- guided 

sin k„ = 0 
n 

nir 

Clamped- clamped 

k 

sin ^ = 0 
2 

2n7T 
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I 



(a) Coordinate system for uniform vibrating beam (b) sign convention for forces and moment 

under constant axial load. acting on arbitrary element. 

Figure 1. - Vibrating beam with arbitrary boundary conditions. 



(a) Coordinate system for cantilevered beam with (b) Forces acting at free end of cantilevered beam 

tip mass under constant axial load. under constant axial load. 


Figure 2. - Cantilevered beam with tip mass under constant axial load. 



(a) Coordinate system for cantilevered beam with 
tip mass under directed axial load. 

Figure 3. - Cantilevered beam with 



(b) Forces acting at free end of cantilevered beam 
under directed axial load. 


tip mass under directed axial load. 
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ii IX 



Frequency parameter, pj 

























Figure 14, - Frequency as function of cantilever beam with tip mass. 


Frequency parameter. 



Figure 13. - Fundamental frequency of cantilevered beam with tip mass 
under constant axial load. 


NASA- Langley, 1975 E-8433 



Figure 16. - Frequency as function of axial load for cantilevered beam 
with tip mass. 



Figure 15. - Fundamental frequency of cantilevered beam with tip mass 
and axial load directed through root. 
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